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We investigate the breakdown of disordered networks under the action of an increasing external — 
mechanical or electrical — force. We perform a mean-field analysis and estimate scaling exponents 
for the approach to the instability. By simulating two-dimensional models of electric breakdown and 
fracture we observe that the breakdown is preceded by avalanche events. The avalanches can be 
described by scaling laws, and the estimated values of the exponents are consistent with those found 
in mean-field theory. The breakdown point is characterized by a discontinuity in the macroscopic 
properties of the material, such as conductivity or elasticity, indicative of a first order transition. 
The scaling laws suggest an analogy with the behavior expected in spinodal nucleation. 

PACS numbers: 62.20.Fe, 62.20.Mk, 64.60.Lx 



I. INTRODUCTION 

The breakdown of solids under external forces is a longstanding problem that has practical and theoretical relevance 
0. The first theoretical approach to fracture mechanics dates back to the twenties with Griffith's theory (2| which 
says that cracks grow or heal, depending on whether the external stress prevail over the resistance at surface of the 
crack. Since the work of Griffith, a great effort has been devoted to experimentally test the validity of the theory and 
to extend it to various crack geometries and boundary conditions The Griffith theory, in spirit, is very similar to 
the classical theory of nucleation in first order phase transitions |^ . In bubble nucleation, a critical droplet will form 
when the change in free energy due to the bulk forces exceeds that of the surface terms. 

The analogy between first order transitions and fracture has been investigated further by numerical model and the- 
oretical calculations. Several authors suggested that the breakdown point in thermally activated fracture is analogous 
to a spinodal point. Spinodal nucleation contrary to classical nucleation, is characterized by scaling properties and 
fractal droplets. Rundle and Klein analyzing a Landau-Ginzburg equation for the growth of a single crack, showed 
that the system obeys scaling laws expected for spinodal nucleation. Selinger et al. (7]|| have studied the problem in 
mean-field theory and by numerical simulation. They conclude that a solid under stress is in a metastable state and 
when the external stress is raised beyond a critical value, corresponding to the spinodal point, the system becomes 
unstable. The nature of the nucleation process in a stressed solid was studied by Golubovic et al. Jj| using Monte 
Carlo simulations. Recently, in the framework of elastic theory, it was shown JlOf that the point of zero external stress 
corresponds to the condensation point in gas- liquid first order transitions. One of the ambitious goals of these and 
other studies is to formulate a statistical thermodynamics of fracture. 

Most of the theoretical studies we have discussed deal with the situation in which fracture is thermally activated 
and quenched disorder is irrelevant. In most realistic situations, however, the solid is not homogeneous and disorder, 
in the form of vacancies or microcracks, strongly affects the nucleation process HI). For example, cracks may start 
from different defects and coalesce [jll), in contrast with the assumptions of Griffith- like theories. There are situations, 
encountered for example in material testing, in which the system is driven by an increasing external stress and the 
time scale of thermal fluctuations is larger than the time scale induced by the driving. In those cases, the system can 
be effectively be considered as being at zero temperature and only quenched disorder is relevant. This is the situation 
we investigate in this paper. It is also worth to emphasize that the breaking process is in most cases irreversible, so 
opening and closing a crack is not like flipping back and forth a spin. 

The understanding of the breakdown of disordered systems has progressed to a large extent with the use of large- 
scale simulations of lattice models jl2| . In these models a conductor is represented by a resistor network and a elastic 
medium by a spring network or other more complex discretizations. The disorder is modeled by random failure 
thresholds or by bond dilution. In this way the model retains the long range nature and the tensorial structure of the 
interactions, which are computed solving coupled linear equations. These models have provided a good description 
of geometrical and topological properties of fracture, leading to the injection in this field of scaling concepts Jl3| . 
Recently, quasi-static lattice models have been used also to study dynamical properties of fracture jl4|-|l7|| . 
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The breakdown of a disordered solid is preceded by intense precursors in the form of avalanches. It has been 
experimentally observed that the response (acoustic emission) to an increasing external stress takes place in bursts 
distributed over a wide range of scale. Examples are found in the fracturing of wood p8[ , cellular glass and 
concrete ^0|, in hydrogen precipitation El), in dislocation motion in ice crystals |22] and in volcanic activity |p3| . 
These phenomena are reminiscent of the Gutenberg-Richter law for earthquakes statistics p3] . It is becoming apparent 
that avalanche response is rather the rule than the exception in driven disordered systems. Other examples range from 
the motion of domain walls in magnets (the Barkhausen effect) |25| and flux lines in superconductors |26|], frictional 
sliding (^?J, to fluid flow in porous media |^8| and the inflation of degased lungs J2{J. Therefore, understanding the 
general physical mechanisms of avalanche dynamics goes well beyond the study of breakdown and fracture. 

In this paper, we numerically study the random fuse model po|"|3^| and a spring network jj4j. The investigation 
of these two models allows us to compare the behavior of a quasi-static scalar model, with a more complex vectorial 
dynamical model. We analyze the scaling close to the breakdown and we find that it is consistent with a mean-field 
analysis. We show that avalanche behavior near the breakdown in disordered systems is analogous to the formation 
of droplets observed close to a spinodal instability in first-order phase transitions. The system is driven by a slowly 
increasing external force through a complex energy landscape; it is not allowed to jump over energy barriers by thermal 
activation. This situation should correspond to the experiment reported in Ref. Raj. The intriguing consequence of 
this analysis is that the behavior of a disordered driven system at zero temperature is similar to what is expected from 
a thermally driven homogeneous system close to a spinodal point. It is tempting to conclude that quenched disorder 
has an effect similar to thermal fluctuations, although a discussion in terms of metastability and nucleation is not 
possible in the first case. We will briefly discuss these analogies in the magnetic context. Finally, it is interesting to 
remark that spinodal nucleation and first-order transitions have also been suggested to play an important role in the 
physics of frictional sliding and earthquakes |35| . Also for these systems thermal disorder is expected to be irrelevant 
with respect to quenched inhomogeneities. 

The paper is organized as follows. In section II, we briefly review spinodal nucleation and we discuss the role of 
quenched disorder. In section III, we present a mean-field analysis of fracture. Section IV discusses the simulations 
for the random fuse model and section V is devoted to molecular dynamics simulations of a spring network. In section 
VI, we discuss the cluster structure of our models and draw analogies with percolation and in section VII, we compare 
the results of our model with experiments and discuss some open questions. 

A short account of a subset of these results has been presented in Ref. |l5| . 



II. SPINODAL NUCLEATION, THERMAL FLUCTUATIONS AND QUENCHED DISORDER 

Nucleation near a spinodal appears to be very different than classical nucleation and the classical theory is expected 
to fail. Droplets appear to be fractal objects and the process of nucleation is due to the coalescence of these droplets, 
rather than the growth of a single one |3q ] . The theoretical description of homogeneous spinodal nucleation is based on 
the Landau-Ginzburg free energy of a spin system in presence of an external magnetic field Q . When the temperature 
is below the critical value, the free energy has the typical two-well structure. In the presence of an external magnetic 
field, one of the wells is depressed with respect to the other, which represents therefore the metastable state. The 
system must cross a free energy barrier to relax into the stable phase. When the external field is increased, this 
nucleation barrier decreases, eventually vanishing at the spinodal. Using this formalism, it has been shown that the 
approach to the spinodal is characterized by scaling laws, analogous to critical phenomena. The magnetization or the 
order parameter <j> scales with the external field H as 

^-^-(Hs-H) 1 / 2 (1) 

where (f> s and H s are respectively the order parameter and the field at the spinodal. This law implies a divergence of 
the quasi-static susceptibility 

X =^~(H S -H)-^, 7 =l/2 (2) 

The fluctuations in the order parameter can be related to suitably defined percolation clusters, whose sizes turn out to 
be power law distributed with an exponent r = 3/2, in mean-field theory. For finite-dimensional short-range models, 
this mean-field picture is expected to fail, since the system will nucleate before reaching the limit of metastability. 
On the other hand mean-field behavior is expected to be valid in presence of long-range interactions, and it has been 
numerically verified in Monte Carlo simulations of long range Ising model |}7j . The limit of stability in thermally 
activated homogeneous fracture is believed to correspond to a spinodal point. One should then be able to observe 
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scaling laws consistent with those found in spinodal nucleation. To our knowledge, such a scaling behavior has not 
yet been observed in numerical simulations. 

In this paper we will not deal with thermally-activated fracture but rather with a disordered driven system. In this 
regard, an interesting analogy can be made with a model recently proposed by Sethna et al. f3§j , |3l| in the context of 
magnetic hysteresis. The model in question is a random- field Ising model (RFIM), driven at zero temperature by an 
increasing uniform magnetic field H. Each spin Sj takes the sign of the local force 

fi = -g^ = J^Sj + hi + H, (3) 

where the sums runs over the neighboring sites, and hi is the random field at site i, which has a Gaussian distribution 
with variance R. When the external field is increased, some local forces change sign and the spins flip along the 
direction of H in avalanches. For low values of the strength of the disorder R < R c , there is a critical value of the field 
H c for which the system undergoes a discrete (first-order) transition involving a finite change in the order parameter 
(the magnetization). 

In mean-field theory, the approach to the instability H Cl is characterized by the same scaling laws and exponents 
of spinodal nucleation, as reported in Eqs. (Q)-(j^). Close to the first-order transition, the avalanche size distribution 
is described by a scaling form 

P(m)~ a -Tf[m(H e -H) K ], (4) 

with r = 3/2 and k = 1. It is worth noting that these scaling exponents coincide with the mean- field exponents 
for the distribution of droplets in homogeneous spinodal nucleation. From these studies it appears that the behavior 
of thermally activated homogeneous spinodal nucleation is similar to the approach to the instability in disordered 
systems driven at zero temperature. However, one should bear in mind that for a given realization of the disorder the 
dynamics is completely deterministic in the second case. Concepts such as metastability and nucleating droplets are 
formally not defined in this context. 



III. MEAN-FIELD THEORY 

In this section we generalize the analysis of Ref. B to derive a simple mean-field theory for fracture. The models 
we will analyze are defined on a two-dimensional lattice. Each bond of the lattice is supposed to obey the equations 
of linear elasticity, until it is stretched beyond a randomly chosen threshold, after which it breaks. In the electric case 
the equations are scalar, each bond satisfies the Ohm law, and the currents are computed numerically by solving the 
Kirchhoff equations with the appropriate boundary condition. 

To illustrate the mean-field theory we will consider for simplicity the random fuse model. To every bond i of the 
lattice we associate a fuse of unit conductivity <7j = 1. An external current / or voltage V is then applied to the 
system by imposing an external voltage V to two opposite edges of the lattice. When the current in the bond exceeds 
a randomly distributed threshold Di the bond becomes an insulator (<x; = 0). The voltage drops (AV), for each bond 
are computed by minimizing the total dissipated energy 

£({a})^J>[(An 2 -A 2 ]- (5) 

i 

The dynamics of the model results from a double minimization process. The voltage drops (AV)i are obtained by a 
global minimization of the energy at fixed er^, while the o~i are then chosen to minimize the local bond energy. The first 
step is equivalent to solving the Kirchhoff equations for the network, while the second step corresponds to breaking 
the bonds for which the current overcomes the threshold. The external current is increased slowly until the lattice is 
no longer conducting. This means that each time a bond is broken, the voltage and the currents are recomputed with 
the new values of the conductivities. 

To derive a mean-field theory, it is useful to recast the dynamics of the model in terms of the externally applied 
current /. We can rewrite the energy of Eq. (||), in full generality, as 

where G({a}) is the total conductivity of the lattice and is a complicated function of the local conductivities. We can 
estimate G({a}) using the effective medium theory [Q, which in our case gives 
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G({<x})=20-1, (7) 
where (j> = J^i ^i/L 2 - We can express the energy as a sum over "spins" interacting with effective random fields hi 



E.uf(/,W) = 5>* = f f u£ , - A 2 ) • (8) 



2 ^ V L2< M 2< ^- 1) 

the value of <p can be computed self-consistently as 



1-I/lLy/ '0(20-1)] 

= P(/n < 0) = 1 - / |0(£>)d-D, (9) 



o 

where p(-D) is the distribution of failure thresholds. The solution of this equation can be expressed in terms of the 
current per unit length / = I/L. We can identify / with the external field and <f> with the order parameter. 

We can show (see Appendix) that under general conditions p in Eq. (||) has a solution for f < f c and, close to f c , 
(f> scales as 

<j>-0c~Uc-f) 1/2 - (10) 

The mean-field theory we have presented is very similar to the fiber-bundle model (FBM) with global load sharing, 
an exactly solvable model for fracture which has been studied extensively |ll],[42]]. In the FBM an external load F, is 
applied to N parallel fibers, and equally shared among the unbroken ones. This means that the force on each fiber is 

fi = F/n (11) 

where n = N<f> is the number of unbroken fibers. A fiber breaks when its force exceeds a quenched random threshold 
D. One can write an equation for the density of unbroken fibers that has the form of Eq. (||), with the upper limit of 
integration replaced by F/(N<fr). The FBM can be obtained as a mean- field theory in the case of site damage, since 
in this case the effective medium conductivity is given by G(cj>) ~ </>. 

We can obtain the mean-field avalanche size distribution from the exact results derived for the FBM jl2j 

P(m) ~ m- T /M/ c - /)"); r=| K = 1, (12) 

where m is the number of bonds that break as function of the current. Eq. (|lj) can also be obtained in the case of 
bond damage using similar arguments. 

The average avalanche size (m) is proportional to the "susceptibility" d<j)/df Q, and therefore diverges at the 
breakdown as 

M~(/c-/r 7 7=1/2- (13) 

The exponents we have introduced satisfy the scaling relation k(2 — r) = 7, which is consistent with the values reported 
in Eq. ( [l2"| ) and Eq. ( |l3| ) . The mean-field analysis indicates that the system is undergoing a first order transition since 
the order parameter has a discontinuity and the conductivity at f c has a jump from G{4> c ) > to zero. The approach 
to this transition is characterized by avalanches of increasing size, diverging at the transition. 

A similar behavior with the same scaling exponents is observed in the mean-field theory of the driven RFIM [j38|j39[ 
for small disorder. In the RFIM, one observes also a second order transition as the width of the disorder is increased. 
A similar transition does not seem to be present in our system, at least not in the mean- field treatment. It is also 
interesting to note that the same scaling laws describe metastable systems close to a spinodal point. The quasistatic 



susceptibility diverges as in Eq. (13) and droplets are distributed according to Eq. (ftz 



IV. SCALING BEFORE BREAKDOWN IN THE RANDOM FUSE MODEL 



An important issue to address at this point is the validity of mean-field results in the case of real low-dimensional 
systems. It is known that scaling does not hold close to the first-order transition for short-ran ged RFIM in dimensions 
d = 2, 3 |38|,[59| . Similarly, spinodal singularities are observed when interactions are long-range |37j . Elastic interactions 
are intrinsically long-range, which could lead to mean-field behavior even for low dimensions, as we will next show 
numerically 
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We simulate the random fuse model (3(| on a tilted square lattice, with periodic boundary conditions in the 
transverse direction. As we discussed in Sec. Ill, the current in each bond is obtained solving Kirchhoff equations, 
which we do numerically using a multigrid [E3 relaxation algorithm with precision e = 10~ 10 . The distribution of 
thresholds is chosen to be uniform in the interval [1 — A, 1 + A] 0|. We impose an external current I through the 
lattice and we increase it at an infinitesimal rate. When a bond fails, we recompute the currents to see if other failures 
occur. The process is continued until a path of broken bonds spans the lattice and no current flows anymore. 

We determine the cluster size distribution n(s,I), which is defined as the number of clusters formed by s neighboring 
broken bonds when the applied current is /. The moments (Mk(I) = J s k n(s, I)ds is the fc-th moment) of n(s,J) 
describe much of the physics associated with the breakdown process. We determine n(s,I) by averaging over the 
various threshold distribution configurations. The first moment M\{I) is the total number of broken bonds due to 
the current / and is therefore proportional to (/). According to our mean-field picture, the average (m) of the quantity 
m = dMi(I) I dl should then diverge close to the breakdown as (I c — I) -1 / 2 . In Fig. |l|, we plot the number of bonds 
m that break for a given value of the current, in a particular realization of the process. We see that the breakdown is 
highly inhomogeneous with avalanches of increasing amplitude. In order to test the scaling, we plot in Fig. || (rn)~ 2 
as a function of the reduced current I j 1 I c , where I c is the average breakdown current, and we see that the graph is 
linear. 

We also measure the distribution of avalanche sizes, integrating over all the values of the current. The mean- field 
analysis predicts that P{m,I) ~ to~ 3 / 2 f(m(I c — I)) which yields an exponent r' = 5/2 when the distribution is 
integrated over the current. In fact we see that our data are consistent with this exponent (Fig. (^)). We have also 
checked that the cut-off of the distribution increases with the system size. A similar result (r ~ 2.7) for smaller lattice 
sizes (L ~ 40) was previously reported by Hansen and Hemmer | [46|] , who also pointed out the similarity with the 
predictions of the FBM. 

V. SCALING IN THE FRACTURE OF A SPRING NETWORK 

The study of fracture in elastic medium is carried out by molecular dynamics (MD) simulation on a lattice model 
with elastic restoring forces. Our model consists of a L x L (L = 20 and 50) square network with central and 
rotationally invariant bond-bending forces. The potential energy of the network is ]34| 

V = \ X/ + \ ( 5e vk?9ij9jk, (14) 

{ij) {ijk) 

where Srij is the change in the length of the spring between the nearest neighbor sites < ij > from its equilibrium 
value (taken to be unity), and SOijk is the change in the angle between the adjacent springs ij and jk from its 
equilibrium value tt/2 which is taken to ensure the square lattice structure of the lattice at equilibrium, giy = 1 if 
the spring ij is present and is otherwise, a and b are the force constants of the central and the bond-bending force 
terms respectively. In terms of arbitrary length and time scales Iq and to, the equations of motion in dimensionless 
variables involve the parameters Ai = at^jra and A2 = bt\jml\, where m is the mass associated with the lattice sites. 
We choose Z to be the lattice spacing (most naturally). The ratio A 2 /Ai = b/al^ is then a characteristic of the system 
under consideration. We choose Ai = 1 and A 2 = 0.1. The small value of A2, much less than the value of Ai, allows 
the fracture to develop without much deformation of the network. We start with all the springs intact so that = 1 
for all neighboring ij's and with each spring we associate a random breaking threshold D^, chosen from a uniform 
distribution D G [0,2]. 

We impose a constant external force F on the sites of the boundary and the system is allowed to evolve dynamically 
using Ver let's algorithm |34j] , 

n(t + At) = 2n(t) -n(t- At) + F^At) 2 . (15) 

Here, Fi is the force (as determined from Eq. |l4|) and fl(t) is the position vector of the site i at time t. The simulation 
involves discrete time t in steps of At. If the simulation runs for n iterations, then the elapsed time is t = nAt 
while the real time elapsed is nt n At. One way to speed up the relaxation process would be to choose a large value 
for At. However, there is an upper limit to this value for the iteration process to remain numerically stable. This 
limit is proportional to the convergence time for the fastest developing components of the stress distribution, which 
is generally very small in disordered systems. We choose Ai = 0.01. In addition, we add to the evolution a small 
viscous force to damp out excessive oscillations. In the course of evolution, if any spring ij stretches beyond its cut-off 
value Dij, the spring snaps irreversibly and gij for that spring is set to zero. We increase the external force F by 
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small steps and at each step we compute the number of broken bonds, which constitute an avalanche. To average 
over disorder, the simulation is repeated for 50 different configurations of threshold values D. 

The fracture in the network takes place in a series of bursts of bond (spring) breaking processes. In such a burst, 
bonds break from different parts of the network and the fracture grows. We keep track of the clusters formed by the 
connected broken bonds which, when spanning the network, causes its macroscopic breakdown. 

We study here the same quantity analyzed in the previous section, namely the susceptibility and the avalanche 
distribution. In Fig. ||, we plot (m)~ 2 as a function of the stress and we find the linearity as expected. The critical 
breaking stress F c can readily be found from the points of intersection. Next, we consider the distribution P(m) of 
the value of m integrated over all the values of a up to F c . We expect the behavior P(m) ~ m~ 5 / 2 in the mean- field 
theory. Fig. compares our simulation results with the mean-field prediction, The data presented in the figure are 
binned, where for different neighboring m-values, the corresponding P(m)-values are combined in one bin, and the 
result is plotted as the arithmetic mean of the two extreme m-values in the bin. 

VI. CLUSTER ANALYSIS: NUCLEATION OR PERCOLATION? 

The avalanche size defined in the previous sections does not represent the geometric structure of the cracks, but 
only counts the number of bonds breaking at each time step. A geometrical characterization of the damage can be 
obtained studying the cluster size defined counting the number of connected broken bonds. 

In the fuse model, the average cluster size S = M2/M1 increases with I. However, by plotting S for different system 
sizes, we observe that the cluster size is not diverging (Fig. ^). To clarify this point, we confirm that S(I C ) does not 
show scaling with the lattice size L. We find similar results for the spring network, where the cluster size distribution 
has an exponential cut-off that does not change with the lattice size. We also study the number of clusters n c = Mo 
as a function of the current and for different system sizes (see Fig. ^|). We observe that n c scales as 

n c = L 2 g(I/L), (16) 

which is expected for a first order transition. We obtain a similar scaling plot for the spring network (see Fig. ^o|) . 

Next we study the behavior of the lattice conductivity in the fuse model. For a given realization of the disorder the 
conductivity has a discrete jump at the breakdown (Fig. (0a). We also plot the conductivity averaged over different 
realizations of the disorder and we observe a smooth curve with a slope at the breakdown that becomes sharper as 
the system size increases (Fig. @b). In the spring network, we calculate the lattice elasticity Y, which shows similar 
behavior as a function of the applied stress (Fig. ||). 

Two principal scenarios have been proposed to explain the scaling behavior of avalanches prior to rupture. The 
first scenario invokes a continuous phase transition with a diverging characteristic length p7| . The various cracks 
inside the lattice should grow until one of them finally rules over the others, becoming the incipient spanning cluster. 
This is precisely what happens in percolation when the occupation probability p is increased toward the percolation 
threshold p c . If this scenario is true for fracture, we would expect the cluster characteristic size to diverge, contrary 
to our results. In the random fuse network a percolation transition is expected only in the limit of infinitely wide 
disorder distributions |4f|, when the strength of the disorder clearly dominates over the interactions. 

The second scenario, in favor of which we presented numerical and theoretical evidences, describes fracture as a 
first-order phase transition close to a spinodal-like instability. The clastic state is considered to be metastable, as 
soon as a non-zero stress is applied. Due to the presence of disorder, the system evolves through a series of metastable 
states towards the final instability. This occurs with the nucleation of cracks growing up to a critical size s c at which 
they coalesce forming the macroscopic crack. Contrary to percolation, in this case there is no incipient spanning 
cluster prior to rupture. 

What explains then the scaling in the avalanche statistics and the susceptibility? We recall that elastic (or electric) 
forces are long-range. When nucleation occurs close to a spinodal instability — which is well defined only for mean-field 
or long-range interactions — one expects a divergent susceptibility Q . This is not naively related to the fluctuation of 
a geometrical quantity such as the crack size, which is not diverging at the spinodal. In order to describe geometrically 
the susceptibility, it is necessary to define the clusters in a peculiar way, considering each site connected with all the 
others within the range of interactions |37f| . These fluctuations are therefore different from those encountered in a 
second-order phase transition. 

The spinodal point is a quite peculiar critical point which, rigorously speaking, exists only in mean-field theory, 
but can be detected when long-range interactions are present. In this respect, no scaling is observed in the avalanche 
distributions when the stress transfer after breaking is local, such as in the local load-sharing fiber bundle model 
studied in Ref. @g§. 
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VII. EXPERIMENTAL COMPARISON AND OPEN PROBLEMS 



The above results clarify the nature of the breakdown process in the presence of quenched disorder. We have found 
that the breakdown is preceded by avalanches distributed as power laws. The scaling exponents are in quantitative 
agreement with the prediction of mean-field calculations. We have discussed that only globally defined quantity such 
as (m(/)) and P(m) display scaling, while locally defined quantities, such as S, do not show any singular behavior. 
For a second-order phase transition, we would expect local quantities to show scaling. For instance, in percolation 
we have S ~ (p ~ Pc) 1 , where p is the concentration of broken bonds. On the other hand, first-order transitions 
usually do not show any precursor and scaling is not observed. An exception to this rule is represented by first-order 
transitions close to a spinodal point, for which some global quantities display scaling [ j37| ; we will argue that this case 
is relevant to the behavior observed before breakdown (49| . 

The observation that mean-field scaling is present in the fracture of two different network models suggests that 
this behavior is rather robust and does not depend on the fine details of the models, such as the tensorial structure 
of the interactions, the dynamics or the boundary conditions. A more stringent test of our conclusions should come 
from the analysis of experimental data. In particular, the experimental setup discussed in Ref. fl8|| resembles some of 
the features of our model. An external pressure is slowly increased until the material (wood or fiberglass) breaks. In 
this process acoustic energy is released in bursts, whose amplitude shows a net increase as the material approaches 
the breakdown point |jq ]. The integrated distribution of burst energies E was found to follow a power law with 
an exponent roughly equal to —2, which must be compared to r = 5/2, if one assumes that m in our paper is 
proportional to E in Ref. Jl||. We see that there is a discrepancy in the results, although this may be due to the 
statistics. We have also tried to analyze the scaling of the average energy released at pressure P, but we could not 
obtain a firm conclusion due to the large statistical uncertainties. Interesting results have been recently obtained in 
three dimensional simulations of fuse networks |5(| . 

In mean-field theory, driven disordered systems behave similarly to their homogeneous, thermally driven, counter- 
parts, if we compare the scaling of avalanches with that of the droplets. This applies to the RFIM which shows 
features similar to those of spinodal nucleation [^|, and to the fracture models we have studied. However, one should 
be careful not to interpret these analogies too strictly, since in driven disordered systems the notions of metastability, 
spinodal point and nucleation are not well defined. In particular, the identification of 4> and / with the order and 
control parameters is justified only in MF theory. In two dimensions, simple homogeneous scaling fails in the presence 



of disorder. The breakdown current I c has a logarithmic size dependence |31 



Ic-T^J, (17) 

which cannot be interpreted in MF theory. A similar dependence is also present in the fraction of broken bond before 
breakdown </> c . While for finite systems the picture we have presented is completely consistent, it is not obvious how 
to perform the L — > oo limit. In order to obtain intensive parameters in this limit one should rescale / and (f> by 
an appropriate logarithmic factor. This can be done implicitly analyzing the data in terms of f/f c , as in Ref. [ jl8| . 
Finally, we note that in most cases the final breakdown starts from existing defects in the material. In our simulations, 
we restricted our attention to the case in which these defects were smaller than the discretization unit. 

In conclusion, we have shown that two different models of breakdown and fracture share the same mean-field 
scaling exponents approaching the rupture point. The behavior observed in these models is analogous to spinodal 
nucleation in thermally driven homogeneous systems. At the breakdown point, the macroscopic quantities (elasticity, 
conductivity) are discontinuous and the characteristic crack size S stays finite in the large L limit. In addition, 
the statistics of global quantities (i.e. the number of broken bonds) display clear mean-field scaling in analogy with 
spinodal nucleation. The direct applications of these ideas to experiments still remains an open question. 
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APPENDIX 



We derive here the scaling law for the approach to the instability in mean-field theory. We define f = I/L and 
h((f>) = y / 4'(2(j) — 1) and we rewrite Eq. as 

4>=1- p(D)dD. (18) 



(19) 



By taking the derivative over / on both sides of the equation, we obtain for the susceptibility 

d<f> = Pif/hm 

df {i- P {f/h(4>))fh>/h(4>) 2 Y 



When the denominator is equal to zero, the system reaches the instability and the susceptibility diverges, which 
defines the critical values <f> c and f c 

(1 - p(f c /h(0 c ))f c ti(<f> c )/h(<t> c f) = 0. (20) 
The Taylor expansion of Eq. (|l^) around (<f> c , f c ) yields 

8<j> = 4> - 4> c = p{fclh{4> c )){f c ti{4>c)mh{4> c ) 2 

-5f/h^ c ))+A5cf> 2 +B8f8cp, (21) 



where Sf = (/ — f c ) and A, B are two constants. The term proportional to 5<f> vanishes because of Eq. (|20|), leaving 
an equation of the form 

8f ~ 8<f> 2 (22) 

which is the scaling relation reported in Eq. (|To|). This relation thus holds for any analytic normalizable distribution 
function, but it is also true in the case of a uniform distribution, as can be easily shown by a direct calculation. 
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FIG. 1. Avalanches in the fuse model. As the current is 
increased the bonds break in avalanches of increasing size until 
the final breakdown occurs. 
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FIG. 2. The average avalanche size in the fuse model scaled 
with the mean-field exponent (7 = 1/2) as a function of I/I c , 
for different values of the system size L. The linearity of the 
plot supports the validity of the mean-field calculations. 
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FIG. 3. The avalanche size distributions in the fuse model 
for two values of system size plotted in log-log scale. A line 
with the mean-field value r' = 5/2 of the exponent is plotted 
for reference. 
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FIG. 4. The average avalanche size in the spring network 
scaled with the mean-field exponent (7 = 1/2) as a function 
of the applied stress F, for two different values of the system 
size L. The linearity of the plot supports the validity of the 
mean-field calculations. 
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FIG. 5. The avalanche size distributions in the spring net- 
work for systems with two values of system size plotted in 
log- log scale. A line with the mean-field value r' = 5/2 of the 
exponent is plotted for reference. 
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FIG. 6. The average cluster size as a function of the current 
for different system sizes. Note that the cluster size does not 
diverge. 
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FIG. 7. (a) The conductivity as a function of the current 
for a single realization of the disorder, for different system 
sizes. (b) The conductivity as function of I/I c averaged over 
different realizations of the disorder. Note that the discrete 
jump, indicative of a first-order transition, is smoothed for 
small system sizes. 
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FIG. 8. The elasticity of the spring network as a function 
of the applied stress. 
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FIG. 9. (a) The number of clusters as a function of the 
current in the fuse model for different system sizes, b) The 
corresponding scaled plot. 
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FIG. 10. The scaled number of clusters as a function of the 
stress in the spring network for different system sizes. 
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